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Abstract 

The technique of Pade Approximants, introduced in a previous work, is applied to 
extended recent data on the distribution of variations of interest rates compiled by the 
Federal Reserve System in the US. It is shown that new power laws and new scaling 
laws emerge for any maturity not only as a function of the Lag but also as a function of 
the average inital rate. This is especially true for the one year maturity where critical 
forms and critical exponents are obtained. This suggests future work in the direction of 
constructing a theory of variations of interest rates at a more "microscopic" level. 
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1 Introduction 



Many models have been constructed to try to describe and/or predict the term 
structure of interest rates. Recently, it has been shown on a specific set of data that 
the empirical distribution of the variation of interest rates is rather well reproduced 
by using the technique of Pade Approximants pQ. This article is a further step 
in this direction. The classical approach is based on Brownian motions of various 
kinds. It leads to exponential tails (or superposition of exponential tails) in the 
distributions. These so called narrow tails are usually badly rejected by the data. 
Levy- flight 0, j3] distributions have also been proposed. In this case the tails 
of the distributions decrease as the inverse of a small power of the variation of the 
interest rates (see below (jl(Jj) ). in fact a power comprised between one and two. The 
predicted tails are really too fat and are also rejected by the data. Sometimes rather 
unjustified cuts have been used to sharpen the Levy tails. 

In a recent paper |T], analyzing numerically the data (Data {1}) on the Ameri- 
can daily spot interest rates compounded by the Federal Reserve System |5] between 
February 15, 1977 and August 4, 1997 (iV events = 5108 corresponding to 5108 open- 
ing days), it was discovered that the fat tails, known to appear in the empirical 
distributions, decrease essentially as a fourth power of the variation of the interest 
rate and are thus naturally amenable to ratios of polynomials. It was shown that 
the empirical distributions, averaged on the initial interest rates, are very well fitted 
with simple Pade Approximants. More recently, articles taking into account the fat, 
but not too fat, tails and power laws have been proposed jHl-jll- 

In this paper, using more complete data (Data {2}) recently available, again from 
the FRS between January 2, 1962 and December 27, 2002 (10230 opening days, see 
(JH1) below), we have been able to study in a more detailed way the Pade fits. At 
a first stage, we confirmed the preceding analysis. The results are summarized in 
section 2. In section 3, we have refined the analysis by dividing the new data in 
suitable subsets, restricting ourselves to intervals in the initial interest rate. This 
has allowed us to confirm and to extend the validity of the preceding results, but 
also and more importantly to suggest new scaling laws. 

These results will be used in a forthcoming article |B] to try to build a "micro- 
scopic model" of interest rates. 
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2 Pade Approximants fits. Update of the Previ- 
ous Results 



In this section, we present briefly the notation and update the results, which have 
been obtained in PQ for the Data {1}, using the full data set (Data {2}) available 
today 

Suppose that, at some initial time to, the interest rate (for some maturity [m]) 
is Iq and, at some later time tf, it is If. The lag and the variation of interest rates 
were defined as follows 

• The lag is the time interval between the initial and final times 

L = t f -t . (1) 

It should be stressed that the lag is expressed in consecutive opening days and 
that non working days were simply discarded j^]. 

• The variation of interest rates is 

v = I f -I . (2) 

Since, in the data, the interest rates are expressed exactly in basis points 
(0.01 percent), the basis point has been chosen as the natural unit. When 
discretized, the variables I and v are thus represented by integers, / and v. 

Let us denote by p^(L, v, Iq) the probability distribution and by p^(L,v, Iq) 
the corresponding discretized distribution. They are the probability that the interest 
rate for maturity [m] has moved from Iq to If during the lag L. Their norms are 

[°° pW(L,v,I )dv = 1, (3) 

J — oo 
oo 

JT rM(L,v,T ) = i. (4) 

v=— oo 

In the FRS Data {1} and {2}, the maturities [m] = 1, 2, 3, 5, 7, 10 and 30 
years are given. All the daily interest rates are expressed exactly in basis points. 
They vary roughly from 1.5 % — 150 basis points to 18 % = 1800 basis points. The 
discrete empirical distributions w(I) of the initial interest rates, for a given maturity 
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and for the two sets, are defined by 
N d (T) 



w(I) 
N d {T) 



N, 



events 



-^events 

number of days the interest rate, for maturity [m] 
was / basis points in the FRS Data {1} or {2} , 
5108 in Data {1} for all [m] , 



(5) 





m 




1 


2 


3 


5 


7 


10 


30 


-^events 


10230 


6606 


10199 


10199 


8323 


10199 


6214 



in Data {2} . 



(6) 



As an example, the FRS empirical distribution w for the maturity [m] = 1 year and 
for the Data {2} is given in Figure (JTJ). Obviously, we have for any such distribution 



1800 



E m = 1 

7=150 



(7) 



The main subject of study in pQ was in fact the discrete average p^l TS , ge (L, v) 
defined by 

1800 

P l ™L gc (L, v) = JT tfM(Io) P [m] (L, v, To) . (8) 

A)=150 

The discretized averages can easily be extracted from the data for lags (in days) such 
that 1 < L < 30 and for variations of interest rates v (expressed in basis points) 
included in the range R v = [—200 < v < 200]. Indeed, outside this range, the data 
is always zero. The discrete normalization becomes 



200 



y y Paverage (-^J V ) 



1 . 



(9) 



d=-200 



It has been first shown empirically, using Data {1}, that the tail of the term 
structure of the interest rates decreases essentially as a fourth power of the variation 
of the interest rates i.e. 



P. 



average 



(L,v) 



K ■ 



when \v\ becomes large 



with a suitable constant k and a power d 

dm 4 . 



(10) 
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This means that the distributions decrease much faster than what Levy flight models 
suggest (namely 1 < d < 2) but much more slowly than what Brownian models 

Paverage ~ « exp (-V 2 / <J 2 ) predict. 

Using detailed fits, it has been shown in pQ that the average distribution (jHJ) can 
be very well (amazingly well) approximated by a Pade Approximant P(0,4), i.e. 
with a constant numerator and a polynomial of fourth degree in the denominator. 
More precisely, for all maturities [m] and for all lags L, the following Pade form of 
the probability distributions, in terms of v, 

Paverage(-^5 V ) = 7~ ~ 7ZZ _.\ ) (12) 

n(l + (q 2 + 2q 2 )v 2 + q%v 4 ) 

fits the Data {2} (and the Data {1} set) rather well. 

Let us enumerate a few of the properties of this form and of the two remaining 
real parameters q\ and q 2 which depend on the lag L and on the maturity [m]. 

1. The continuous distribution corresponding to (JT2J) is, for qi positive, a purely 
positive real function. In fact it is constructed as the modulus square of the 
complex function / 



f[m\ it \ /io\ 

/averaged, ^ (1 + {q lV + q^) > ^ 

P%LjL,v) = | f™ age (L,v) | 2 . (14) 

2. The normalization (J3j) of the distribution ()14|) is automatic by construction. 

3. The distribution depends only on v 2 , hence it is symmetrical under the in- 
terchange {v <-» —v}. This symmetry was checked empirically to hold to a 
very precise degree. In other words, violations of this symmetry are zero well 
within one standard deviation. 

4. The tail of the distribution (|T2|) is of the form v ~ 4 . The Bayesian Hill estima- 
tors PH] of the empirical tails favor an approximate decrease of this kind. 

5. An estimate of q\ is given by 

9i = vrpH rage (L,0) (15) 

and of q 2 by 

q 2 = -a' 1 (16) 
where a stands for the variance of the data distribution. 
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6. The parameters q\ and q 2 have very simple regular forms in term of the lag 
L when L varies from one day to about a month and for any maturity [to]. 
Essentially 



where, for convenience, we have chosen one arbitrary reference value L r = 
15 days. The /ij, which we call scaling parameters, and the less essential 
multiplication constants $ depend on [to] but not on L. The Table (JTJ) (resp. 
Table (0)) shows the /ij and the ratio ^2/^1 that we have obtained for the 7 
investigated maturities using the Data {2} (resp. Data {1}) set. 

In other words, as illustrated by the Figures (J2J) and (jUJ) for [to] = 2 years, as 
an example, in a ln(gj) versus fn(L) plot the data points align along a straight 
line with slope — /ij. The plots for the other [to]'s follow analogous straight 
lines. 

7. Moreover, one may note that one critical exponent is approximatively related 
to the other exponent by (see Tables (JTJ) (J2J)) 



8. It should be stressed that this scaling law (fTTj) seems to be an important 
discovery [I] which had escaped attention. 

3 Pade Approximants fits. Initial Interest Rate 
Dependence 

More recently, we have studied the extended set of data (Data {2}) and performed 
a more refined analysis. 

Indeed, the new data having often a higher statistics, meaningful subsets of 
data can be defined with initial interest rates limited to intervals. We have found 
that subsets containing about two thousand events are large enough to lead to a 
sufficiently precise determination of the parameters. 

A detailed numerical analysis of the distributions restricting the data to an al- 
most fixed initial interest rate Iq has convinced us that the same Pade behavior 
prevails. An approximate determination of the dependence of the Pade parameters 
qi on initial interest rate I follows. It suggests that new scaling laws are at work. 

Let us summarize our results 




1 = 



1,2 



(17) 



fj, 2 ~ 1.8 /ii , 



(18) 
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1. We have first divided the Data {2} set in subsets where Iq is restricted to 
regions (overlapping or not) 

A = [/omin) ] • (19) 

The total probability u A to be inside A and the total numbers of events 
inside A are given by 

max 

= > (2°) 

N A = u A N events . (21) 

In the present analysis, the regions have be chosen in such a way that the 
A^a's are around two thousands. If the number of events is chosen to be 
smaller than about two thousand, the accuracy in the determination of the 
parameters becomes hazardous. On the other hand, there seems to be no need 
to select larger regions with more events in them. The events in Data {2} 
can thus, for example for [m] = 1 year, be separated in five non overlapping 
intervals covering the full Iq range. 

2. Let us define Iq as the empirical average of / in a given region A 

_ w(T) T 

io = 7 ~ J ° min . (22) 

LU A 



3. The distributions, for the events within a region A (see (j2T|) for a precise 
definition), are then supposed, in first approximation, to be represented by 
normalized (jBJ) Pade forms 

p{L, v, I ) = ^ . (23) 

7r(l + (qf + 2q 2 )v 2 + q$v A ) 



4. The parameters % appearing in Eq. (j23j) (for a given maturity [m]) are thus 
evaluated as functions of L and Iq. They are obtained from the data by 
minimizing a x 2 function. In about all instances the fits are very good as 
exemplified by the \ 2 values. 

5. In Figure (0), as an example using the Data {2} with L = 1 and [m] = 1 year, 
two thousand points in the region A =[420 basis points, 553 basis points] 
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(see f)19|l ) have been selected and the empirical points are plotted versus the 
best Pade fit. The average initial interest rate is about J = 502 basis 
points. In this case, the values of evaluated in their respective units are 
gi = 38.6 ± 1.3(%)" 1 and g 2 = -308 ± 19(%)" 2 . The obtained x 2 is 61 with 
83 degrees of freedom. In almost all L, [m] cases an eye check of the curves is 
as satisfactory and the computed x 2 are as good. 

6. At this point, an important remark is worth making. It should be stressed 
that the parameters g, obtained this way are not really those at I = I itself 
but a complicated average of the approximate g^ in the corresponding A range. 
Indeed, it is obvious that the superposition of two or more Pade forms is, at 
best, only approximatively of a Pade form. Since this is an important point, 
let us explain it more precisely. If p(L, v, I) is known for any J, then the 
average probability p(L, v, I ) (J23|) on A is defined by 

_ Y%T ti(T)p(L,v,T) 
p(L, v, Jo) = 7 " Jomin • (24) 

If it then were true that p(L, v, I) is exactly of the Pade form (which is anyhow 
not a correct statement but, at best, an approximate one), the A-average 
probability p(L, v, Iq) would not exactly be of a Pade form. 

7. After the fits have been performed, we find that the new parameters qi depend 
effectively and importantly on I . 

8. In a first approach, we have chosen to minimize the natural x 2 leaving qi as 
free parameter and by determining q 2 from the experimental variance using 
the relation (|T6|) . The interested reader can find the precise definition of the 
X 2 in the reference pQ, in and after their Eq.(22). 

9. We defer to the next section a detailed discussion of the parameters qi, obtained 
from the data, as functions of the lag L and of the average interest rate io- 

4 Pade parameters as a function of the lag, of the 
maturity and, of the initial interest rate 

In this section, we summarize the scaling law and the I dependence that are sug- 
gested by the empirical data. 
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4.1 General considerations and scaling laws 

The first scaling law and the Iq dependence can be summarized as follows. For each 
maturity, the empirical qi are given by the approximate forms (compare to ([17))) 



It) ?< -t • ( 25 ) 




where the Ti are functions depending of the initial interest rate. For convenience, 
we have chosen two arbitrary reference values 

L r = 15 days 

I T = 600 basis points . (26) 

In (I25|) . the lag L is in days, the average initial rate Jo is in basis points, the 
parameters qi and f3\ are in inverse basis points, the parameters q2 and /?2 are in 
inverse basis points squared while the scaling parameters /ij are pure numbers. 

The scaling law (|25|) says that, for every maturity [m], each of the parameters 
qi is the product of two functions. One function depends on the lag L only and has 
the form of a power law as (JT7J). The second function, JF^, depends on I only. 

The best values for the parameters /ij are given in Table ^ as a function of the 
maturity. The form of the Ti can be found in the fourteen plots given in Figures (0)- 
(jSJ). In Figure (jSJ), for [m] = 1 year and L = 1 day, the ln(g x ) is plotted versus Io/I r . 
In Figure (0), for [m] = 1 year and L = 1 day, the ln(g 2 ) is plotted versus Io/I r - In 
Figure (J7|) there are six plots where for L — 1, ln(gi) is given versus Io/I r for [m] = 
2, 3, 5, 7, 10, 30 years. Finally, in Figure (jHJ) there are six plots where for L — 1, 
ln(g 2 ) is given versus I§jl r for [m] = 2, 3, 5, 7, 10, 30 years. 

We have shown that the assumption in (|2"K|) that the qi factorize (separation of 
variables) is favoured by the data. Hence, for a given maturity, the distributions of 
ln(gj) as a function of Iq for different L values just differ by a constant shift. This 
shift on the ln(gj) axis is equal to 

Hqi)-\n(ql)=nML'/L) . (27) 

To illustrate this, the Figures (JH1) and (fTU|) show the distributions of ln(gj), for 
[m] = 5 years as a function of Iq for 3 different L values (1, 10 and 30 days). To make 
the comparison easier, the L=10 days and the L=30 days distributions have been 
moved upwards by their respective shift term (/ijln(10) and ^ ln(30) respectively). 
For the other maturity, one obtains similar plots. 

It has to be emphasized that this variable separation in the approximate forms 
may be a hint of a scaling law which is at work. 
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4.2 The scaling laws for the [m] = 1 year case 



For the [m] > 1 year data, no particular simple shape for the function Ti can be 
extracted. In the [m] = 1 year data, the functions T have a very simple shape. The 
ln(gj) distributions as function of Iq (see Figures ©(jUJ) for L = 1 day) can be fitted 
with straight lines and the following extremely simple approximate forms hold 




ln(/?i) = 3.14 ±0.02, //i = 0.657 ±0.005, v x = 1.32 ± 0.02 (29) 
m(/3 2 ) = 4.51 ± 0.02, /i 2 = 1.197 ± 0.005, u 2 = 2.19 ± 0.02 . (30) 

Hence, for the particular case of [m] = 1 year, we obtain an interplay of two simple 
scaling laws with critical exponents \ii and V{. 



5 Conclusion 

In this paper, we have extended an earlier work aimed at approximating the term 
structure of interest rates by ratios of polynomials called Pade Approximants. The 
form which is used is a purely positive Pade -P[0,4] with a constant numerator and 
a fourth degree denominator in the variation of interest rates. 

We have shown that for any maturity, for any lag and for initial interest rates 
restricted to regions, the form holds to a very good degree of approximation, better 
than what should have been hoped a priori. 

The empirical parameters are represented by functions which point toward the 
existence of scaling laws and scaling exponents at a more "microscopic level" . These 
scaling laws were already discovered and discussed PQ in terms the lag variable. In 
this paper, using the extended FRS data [5] now available, we have shown that 
scaling laws seem also to be at work not only in the lag, as was shown previously, 
but also in the initial interest rate. This is true especially for [m] = 1 year maturity 
where simple critical forms ()28|) appear both in the lag L and in the average interest 
rate Iq and critical exponents can be extracted. 

In a forthcoming paper, we will try to use this result and the scaling laws which 
are suggested by the data to build a theory of the term structure of interest rates 
at a "microscopic level". 
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Figure Caption 



Figure 1. 

The FRS empirical weight distribution w of the interest rates (see (JJJJ) for the 
maturity [m] = 1 year and for the Data set {2}. 

Figure 2. 

Distribution of ln(gi) as a function of ln(L/15) for [m] = 2 years using Data 
set {2}. 

Figure 3. 

Distribution of ln(g 2 ) as a function of ln(L/15) for [m] = 2 years using Data 
set {2}. 

Figure 4. 

One example of a fit to the Data set {2} for [m] = 1 year and L = 1 day 
by a Pade Approximant (continuous curve). The chosen subset (see (JUJ)) 
corresponds to 2000 points in the region Jomin = 420 and / max = 553. The I 
is computed to be 502 basis points. The empirical points are plotted versus 
the best Pade fit. In this case, the values of evaluated in their respective 
units are q x = 38.6 ± 1.3(%) _1 and q 2 = -308 ± 19(%y 2 . The obtained \ 2 is 
61 with 83 degrees of freedom. 

Figure 5. 

The parameter ln(q'i) as a function of Jo/600 for [m] = 1 year and L = 1 
day. These points have been obtained by minimizing the x 2 as a function of 
parameter qi . Each point corresponds to a subset of data containing close 
to two thousand points with average J value equal to I . The subsets are 
overlapping so that the points are not independent. The unit used for q\ 
is For Iq it is the basis point. The straight line is a linear fit. Its 

equation is ln(^) = ln(/5i) - i/i(7o/600) - /iiln(l/15) with v x = 1.32 ± 0.02 
and 1h(#l) = 3.14 ±0.02. 

Figure 6. 
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The parameter ln(g 2 ) as a function of Jo/600 for [m] = 1 year and L = 1 
day. These points have been obtained using the relation defined in Each 
point corresponds to a subset of data containing close to two thousand points 
with average Jo value equal to Jo. The subsets are overlapping so that the 
points are not independent. The error on ln(g 2 ) is computed by propagating 
the error on the estimate of the variance. The unit used for q 2 is (%)~ 2 . For 
Jo it is the basis point. The straight line is a linear fit. Its equation is ln(g 2 ) = 
ln(/? 2 )-i/ 2 (7o/600)-/i 2 ln(l/15) with u 2 = 2.19±0.02 and ln(/3 2 ) = 4.51±0.02. 

Figure 7. 

The parameter In ((ft) as a function of J /600 for [m] = 2, 3, 5, 7, 10, 30 
years and L = 1 day. These points have been obtained by minimizing the 
X 2 as a function of parameter q±. Each point corresponds to a subset of data 
containing close to two thousand points with average Jo value equal to Iq. The 
subsets are overlapping so that the points are not independent. The unit used 
for qi is (%) _1 . For J it is the basis point. 

Figure 8. 

The parameter ln(g 2 ) as a function of J /600 for [m] = 2, 3, 5, 7, 10, 30 years 
and L = 1 day. These points have been obtained using the relation (|T6*j) . Each 
point corresponds to a subset of data containing close to two thousand points 
with average Iq value equal to Iq. The subsets are overlapping so that the 
points are not independent. The unit used for g 2 is (%)~ 2 . For Iq it is the 
basis point. The error on ln(g 2 ) is computed by propagating the error on the 
estimate of the variance. 

Figure 9. 

Distribution of ln((ft), for [m] = 5 years as a function of Jo/600 for 3 different 
L values (1,10 and 30 days) where the L= 10 and 30 days distributions have 
been shifted on the hi^) axis by the correcting term (|2*7j) . 

Figure 10. 

Distribution of ln(g 2 ), for [m] = 5 years as a function of Jo/600 for 3 different 
L values (1,10 and 30 days) where the L= 10 and 30 days distributions have 
been shifted on the hi(g 2 ) axis by the correcting term (|27jl. 
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Table Caption 

Table 1. 

The values of the best fits to the parameters \i{ (see f)25|l ) using Data set {2}. 
These parameters are given as a function of the maturity [m] . The standard 
errors correspond to a unit deviation in the y 2 . 

Table 2. 

The values of the best fits to the parameters \ii (see (J2SJ)) using Data set {1}. 
These parameters are given as a function of the maturity [m] . The standard 
errors correspond to a unit deviation in the x 2 . 
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Figure (|T]) 



weight distribution for [m] = 1 year 




I ^P) 

Figure (0) 



ln(qj) versus ln(L) ([m] = 2 years) 
4 I 

A 

3.5 - 
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Figure (J3J) 



ln(q 2 ) versus ln(L) ([m] = 2 years) 




-3 -2.5 -2 -1.5 -1 -0.5 0.5 1 

ln(L/15days) 

Figure flU) 

A fit for [m] = 1 year, L = 1 day and average L= 502 bp 



300 - 




15 



Figure fl^D 

ln(qj) versus L (L= 1 day & [m] = 1 year) 
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Figure © 

ln(q 2 ) versus I Q (L= 1 day & [m] = 1 year) 
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Figure © 

ln(qj) versus I ([m] = 5 years) 
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Figure (TIT)!) 

ln(q 2 ) versus I„ ([m] = 5 years) 
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Table flU 



Data set {2} 



maturity (yrs) 






- 1 ) 


M2 




- 2 ) 




1 ( 




1 


0.657 


± 


0.005 


1.197 


± 


0.005 


li 


\2 


± 


0.02 


2 


0.627 


± 


0.006 


1.178 


± 


0.007 


li 


18 


± 


0.03 


3 


0.653 


± 


0.005 


1.164 


± 


0.006 


1.78 


± 


0.02 


5 


0.638 


± 


0.005 


1.139 


± 


0.005 


1.79 


± 


0.02 


7 


0.626 


± 


0.005 


1.111 


± 


0.005 


1.78 


± 


0.02 


10 


0.636 


± 


0.005 


1.100 


± 


0.005 


1.73 


± 


0.02 


30 


0.578 


± 


0.006 


1.068 


± 


0.006 


li 




± 


0.02 








Table 
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Data set {1} 














maturity (yrs) 
















1 ( 
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0.607 


± 


0.006 


1.209 


± 


0.009 


1.99 


± 


0.03 


2 


0.644 


± 


0.004 


1.200 


± 


0.007 


li 


56 


± 


0.03 


3 


0.652 


± 


0.004 


1.172 


± 


0.007 


li 


10 


± 


0.02 


5 


0.622 


± 


0.006 


1.139 


± 


0.005 


li 


13 


± 


0.02 


7 


0.606 


± 


0.006 


1.109 


± 


0.004 


li 


13 


± 


0.02 


10 


0.597 


± 


0.005 


1.093 


± 


0.005 


li 


13 


± 


0.02 


30 


0.588 


± 


0.005 


1.074 


± 


0.004 


li 


13 


± 


0.02 
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Emergence of Power Laws and Scaling Laws 
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Abstract 

The technique of Pade Approximants, introduced in a previous work, is applied to 
extended recent data on the distribution of variations of interest rates compiled by the 
Federal Reserve System in the US. It is shown that new power laws and new scaling 
laws emerge for any maturity not only as a function of the Lag but also as a function of 
the average inital rate. This is especially true for the one year maturity where critical 
forms and critical exponents are obtained. This suggests future work in the direction of 
constructing a theory of variations of interest rates at a more "microscopic" level. 
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1 Introduction 



Many models have been constructed to try to describe and/or predict the term 
structure of interest rates. Recently, it has been shown on a specific set of data that 
the empirical distribution of the variation of interest rates is rather well reproduced 
by using the technique of Pade Approximants [1]. This article is a further step 
in this direction. The classical approach is based on Brownian motions of various 
kinds. It leads to exponential tails (or superposition of exponential tails) in the 
distributions. These so called narrow tails are usually badly rejected by the data. 
Levy- flight [2], [3], [4] distributions have also been proposed. In this case the tails 
of the distributions decrease as the inverse of a small power of the variation of the 
interest rates (see below (10)), in fact a power comprised between one and two. The 
predicted tails are really too fat and are also rejected by the data. Sometimes rather 
unjustified cuts have been used to sharpen the Levy tails. 

In a recent paper [1], analyzing numerically the data (Data {1}) on the Ameri- 
can daily spot interest rates compounded by the Federal Reserve System [5] between 
February 15, 1977 and August 4, 1997 (iV cven ts = 5108 corresponding to 5108 open- 
ing days), it was discovered that the fat tails, known to appear in the empirical 
distributions, decrease essentially as a fourth power of the variation of the interest 
rate and are thus naturally amenable to ratios of polynomials. It was shown that 
the empirical distributions, averaged on the initial interest rates, are very well fitted 
with simple Pade Approximants. More recently, articles taking into account the fat, 
but not too fat, tails and power laws have been proposed [6] -[7]. 

In this paper, using more complete data (Data {2}) recently available, again from 
the FRS between January 2, 1962 and December 27, 2002 (10230 opening days, see 
(6) below), we have been able to study in a more detailed way the Pade fits. At 
a first stage, we confirmed the preceding analysis. The results are summarized in 
section 2. In section 3, we have refined the analysis by dividing the new data in 
suitable subsets, restricting ourselves to intervals in the initial interest rate. This 
has allowed us to confirm and to extend the validity of the preceding results, but 
also and more importantly to suggest new scaling laws. 

These results will be used in a forthcoming article [8] to try to build a "micro- 
scopic model" of interest rates. 
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2 Pade Approximants fits. Update of the Previ- 
ous Results 



In this section, we present briefly the notation and update the results, which have 
been obtained in [1] for the Data {1}, using the full data set (Data {2}) available 
today [5]. 

Suppose that, at some initial time t , the interest rate (for some maturity [m]) 
is Iq and, at some later time tf, it is //. The lag and the variation of interest rates 
were defined as follows 

• The lag is the time interval between the initial and final times 

L = t f -t . (1) 

It should be stressed that the lag is expressed in consecutive opening days and 
that non working days were simply discarded [9]. 

• The variation of interest rates is 

v = I f -I . (2) 

Since, in the data, the interest rates are expressed exactly in basis points 
(0.01 percent), the basis point has been chosen as the natural unit. When 
discretized, the variables I and v are thus represented by integers, / and v. 

Let us denote by p^(L, v,I ) the probability distribution and by p' m ](L, v, I ) 
the corresponding discretized distribution. They are the probability that the interest 
rate for maturity [m] has moved from I to If during the lag L. Their norms are 

r pW(L,v,I )dv = 1, (3) 

J — oo 

oo 

jr p [m] (L,v,r ) = i . (4) 

v=— oo 

In the FRS Data {1} and {2}, the maturities [m] = 1, 2, 3, 5, 7, 10 and 30 
years are given. All the daily interest rates are expressed exactly in basis points. 
They vary roughly from 1.5 % — 150 basis points to 18 % = 1800 basis points. The 
discrete empirical distributions w(I) of the initial interest rates, for a given maturity 
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and for the two sets, are defined by 
N d (T) 



w(I) 
N d (T) 



N, 



events 



-^events 

number of days the interest rate, for maturity [to], 
was / basis points in the FRS Data {1} or {2} , 
5108 in Data {1} for all [to] , 



(5) 





m 




1 


2 


3 


5 


7 


10 


30 


-^events 


10230 


6606 


10199 


10199 


8323 


10199 


6214 



in Data {2} . 



(6) 



As an example, the FRS empirical distribution w for the maturity [to] = 1 year and 
for the Data {2} is given in Figure (1). Obviously, we have for any such distribution 



1800 



E m = i • 

7=150 



(7) 



The main subject of study in [1] was in fact the discrete average p^l Tase (L, ^) 
defined by 

1800 

Patage(A v) = ^ W^(I ) §^(L, V, I ) . (8) 
A)=150 

The discretized averages can easily be extracted from the data for lags (in days) such 
that 1 < L < 30 and for variations of interest rates v (expressed in basis points) 
included in the range R v = [—200 < v < 200]. Indeed, outside this range, the data 
is always zero. The discrete normalization becomes 



200 



(9) 



v=~200 



It has been first shown empirically, using Data {1}, that the tail of the term 
structure of the interest rates decreases essentially as a fourth power of the variation 
of the interest rates i.e. 



pivLget-M) ~ K when M becomes large 



with a suitable constant k and a power d 

d w 4 . 



(10) 
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This means that the distributions decrease much faster than what Levy flight models 
suggest (namely 1 < d < 2) but much more slowly than what Brownian models 
^average ~ « exp (-v 2 / a 2 ) predict. 

Using detailed fits, it has been shown in [1] that the average distribution (8) can 
be very well (amazingly well) approximated by a Pade Approximant P(0,4), i.e. 
with a constant numerator and a polynomial of fourth degree in the denominator. 
More precisely, for all maturities [to] and for all lags L, the following Pade form of 
the probability distributions, in terms of v, 

Paveragc(-^> v ) = 7~ ~ ~ 9 _,\ i (12) 

tt(1 + (q 2 + 2q 2 )v 2 + qiv 4 ) 

fits the Data {2} (and the Data {1} set) rather well. 

Let us enumerate a few of the properties of this form and of the two remaining 
real parameters q± and q 2 which depend on the lag L and on the maturity [to]. 

1. The continuous distribution corresponding to (12) is, for qi positive, a purely 
positive real function. In fact it is constructed as the modulus square of the 
complex function / 



/IveLgeOM) " ^ (1 + £ v + g2v2 y (13) 

P [ ™L se (L,v) = | f [ ™L ge (L,v) | 2 . (14) 

2. The normalization (3) of the distribution (14) is automatic by construction. 

3. The distribution depends only on v 2 , hence it is symmetrical under the in- 
terchange {v <-> — v}. This symmetry was checked empirically to hold to a 
very precise degree. In other words, violations of this symmetry are zero well 
within one standard deviation. 

4. The tail of the distribution (12) is of the form v ~ 4 . The Bayesian Hill estima- 
tors [10] of the empirical tails favor an approximate decrease of this kind. 

5. An estimate of qi is given by 

«i = ir$6L.(M) (15) 

and of q 2 by 

q 2 = -a' 1 (16) 
where a stands for the variance of the data distribution. 
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6. The parameters q\ and qi have very simple regular forms in term of the lag 
L when L varies from one day to about a month and for any maturity [to]. 
Essentially 



where, for convenience, we have chosen one arbitrary reference value L r = 
15 days. The //j, which we call scaling parameters, and the less essential 
multiplication constants depend on [to] but not on L. The Table (1) (resp. 
Table (2)) shows the Hi and the ratio H2/ Hi that we have obtained for the 7 
investigated maturities using the Data {2} (resp. Data {1}) set. 

In other words, as illustrated by the Figures (2) and (3) for [to] = 2 years, as 
an example, in a ln(%) versus ln(L) plot the data points align along a straight 
line with slope —Hi- The plots for the other [to]'s follow analogous straight 
lines. 

7. Moreover, one may note that one critical exponent is approximatively related 
to the other exponent by (see Tables (1) (2)) 



8. It should be stressed that this scaling law (17) seems to be an important 
discovery [1] which had escaped attention. 

3 Pade Approximants fits. Initial Interest Rate 
Dependence 

More recently, we have studied the extended set of data (Data {2}) and performed 
a more refined analysis. 

Indeed, the new data having often a higher statistics, meaningful subsets of 
data can be defined with initial interest rates limited to intervals. We have found 
that subsets containing about two thousand events are large enough to lead to a 
sufficiently precise determination of the parameters. 

A detailed numerical analysis of the distributions restricting the data to an al- 
most fixed initial interest rate Jo has convinced us that the same Pade behavior 
prevails. An approximate determination of the dependence of the Pade parameters 
qi on initial interest rate Jo follows. It suggests that new scaling laws are at work. 

Let us summarize our results 




i = 1,2 



(17) 



H2 ~ 1-8 Hi , 



(18) 
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1. We have first divided the Data {2} set in subsets where I is restricted to 
regions (overlapping or not) 

A = [-^Ominj ] • (19) 

The total probability uo A to be inside A and the total numbers of events N A 
inside A are given by 

max 

= > (2°) 

N A = u A N CVCQts . (21) 

In the present analysis, the regions have be chosen in such a way that the 
JVa's are around two thousands. If the number of events is chosen to be 
smaller than about two thousand, the accuracy in the determination of the 
parameters becomes hazardous. On the other hand, there seems to be no need 
to select larger regions with more events in them. The events in Data {2} 
can thus, for example for [to] = 1 year, be separated in five non overlapping 
intervals covering the full Iq range. 

2. Let us define I as the empirical average of / in a given region A 

£%s x w(T) T 
I Q = I - Iomin . (22) 



The distributions, for the events within a region A (see (24) for a precise 
definition), are then supposed, in first approximation, to be represented by 
normalized (3) Pade forms 

tt(1 + [qf + 2q 2 )v 2 + qiv 4 ) 



4. The parameters qi appearing in Eq.(23) (for a given maturity [to]) are thus 
evaluated as functions of L and I . They are obtained from the data by 
minimizing a x 2 function. In about all instances the fits are very good as 
exemplified by the x 2 values. 

5. In Figure (4), as an example using the Data {2} with L — 1 and [to] = 1 year, 
two thousand points in the region A =[420 basis points, 553 basis points] 
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(see (19)) have been selected and the empirical points are plotted versus the 
best Pade fit. The average initial interest rate (22) is about J = 502 basis 
points. In this case, the values of g« evaluated in their respective units are 
q 1 = 38.6 ± 1.3(%) _1 and q 2 = -308 ± 19(%)" 2 . The obtained x 2 is 61 with 
83 degrees of freedom. In almost all L, [m] cases an eye check of the curves is 
as satisfactory and the computed x 2 are as good. 

6. At this point, an important remark is worth making. It should be stressed 
that the parameters qi obtained this way are not really those at Io = Io itself 
but a complicated average of the approximate qi in the corresponding A range. 
Indeed, it is obvious that the superposition of two or more Pade forms is, at 
best, only approximatively of a Pade form. Since this is an important point, 
let us explain it more precisely. If p(L, v, I) is known for any /, then the 
average probability p(L, v, Io) (23) on A is defined by 

_ Ef 0m f w(l)p(L,v,T) 
p(L, v, J ) = . (24) 

If it then were true that p(L, v, I) is exactly of the Pade form (which is anyhow 
not a correct statement but, at best, an approximate one), the A-average 
probability p(L, v, Io) would not exactly be of a Pade form. 

7. After the fits have been performed, we find that the new parameters q t depend 
effectively and importantly on J - 

8. In a first approach, we have chosen to minimize the natural x 2 leaving qi as 
free parameter and by determining q 2 from the experimental variance using 
the relation (16). The interested reader can find the precise definition of the 
X 2 in the reference [1], in and after their Eq.(22). 

9. We defer to the next section a detailed discussion of the parameters q iy obtained 
from the data, as functions of the lag L and of the average interest rate I - 

4 Pade parameters as a function of the lag, of the 
maturity and, of the initial interest rate 

In this section, we summarize the scaling law and the Io dependence that are sug- 
gested by the empirical data. 
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4.1 General considerations and scaling laws 

The first scaling law and the Iq dependence can be summarized as follows. For each 
maturity, the empirical are given by the approximate forms (compare to (17)) 

(25) 

where the Ti are functions depending of the initial interest rate. For convenience, 
we have chosen two arbitrary reference values 

L r = 15 days 

I T = 600 basis points . (26) 

In (25), the lag L is in days, the average initial rate I is in basis points, the 
parameters q\ and l3\ are in inverse basis points, the parameters q2 and i3i are in 
inverse basis points squared while the scaling parameters fii are pure numbers. 

The scaling law (25) says that, for every maturity [to], each of the parameters 
qi is the product of two functions. One function depends on the lag L only and has 
the form of a power law as (17). The second function, T, depends on I only. 

The best values for the parameters fii are given in Table 1 as a function of the 
maturity. The form of the Ti can be found in the fourteen plots given in Figures (5)- 
(8). In Figure (5), for [to] = 1 year and L — 1 day, the ln(gi) is plotted versus Io/I r . 
In Figure (6), for [to] = 1 year and L — 1 day, the ln(g 2 ) is plotted versus Io/I r . In 
Figure (7) there are six plots where for L = 1, ln(gi) is given versus Io/I r for [to] = 
2, 3, 5, 7, 10, 30 years. Finally, in Figure (8) there are six plots where for L — 1, 
In ((72) is given versus Iq/I t for [to] = 2, 3, 5, 7, 10, 30 years. 

We have shown that the assumption in (25) that the % factorize (separation of 
variables) is favoured by the data. Hence, for a given maturity, the distributions of 
ln(gj) as a function of I for different L values just differ by a constant shift. This 
shift on the ln(%) axis is equal to 

hi(Qi)-H<ti=»ML'/L) . (27) 

To illustrate this, the Figures (9) and (10) show the distributions of ln(^), for 
[to] = 5 years as a function of J for 3 different L values (1, 10 and 30 days). To make 
the comparison easier, the L=10 days and the L=30 days distributions have been 
moved upwards by their respective shift term (/Xjln(lO) and /Xjln(30) respectively). 
For the other maturity, one obtains similar plots. 

It has to be emphasized that this variable separation in the approximate forms 
may be a hint of a scaling law which is at work. 
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4.2 The scaling laws for the [to] = 1 year case 



For the [m] > 1 year data, no particular simple shape for the function Ti can be 
extracted. In the [m] = 1 year data, the functions Ti have a very simple shape. The 
ln(gj) distributions as function of I (see Figures (5) (6) for L = 1 day) can be fitted 
with straight lines and the following extremely simple approximate forms hold 




ln(/?i) = 3.14 ±0.02, ix x = 0.657 ±0.005, ^ = 1.32 ± 0.02 (29) 
ln(/3 2 ) = 4.51 ± 0.02, /i 2 = 1.197 ± 0.005, u 2 = 2.19 ± 0.02 . (30) 

Hence, for the particular case of [m] = 1 year, we obtain an interplay of two simple 
scaling laws with critical exponents /ij and z/j. 



5 Conclusion 

In this paper, we have extended an earlier work aimed at approximating the term 
structure of interest rates by ratios of polynomials called Pade Approximants. The 
form which is used is a purely positive Pade P[0,4] with a constant numerator and 
a fourth degree denominator in the variation of interest rates. 

We have shown that for any maturity, for any lag and for initial interest rates 
restricted to regions, the form holds to a very good degree of approximation, better 
than what should have been hoped a priori. 

The empirical parameters are represented by functions which point toward the 
existence of scaling laws and scaling exponents at a more "microscopic level" . These 
scaling laws were already discovered and discussed [1] in terms the lag variable. In 
this paper, using the extended FRS data [5] now available, we have shown that 
scaling laws seem also to be at work not only in the lag, as was shown previously, 
but also in the initial interest rate. This is true especially for [m] = 1 year maturity 
where simple critical forms (28) appear both in the lag L and in the average interest 
rate Iq and critical exponents can be extracted. 

In a forthcoming paper, we will try to use this result and the scaling laws which 
are suggested by the data to build a theory of the term structure of interest rates 
at a "microscopic level". 
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Figure Caption 



Figure 1. 

The FRS empirical weight distribution w of the interest rates (see (5)) for the 
maturity [to] = 1 year and for the Data set {2}. 

Figure 2. 

Distribution of hi(qi) as a function of ln(L/15) for [to] = 2 years using Data 
set {2}. 

Figure 3. 

Distribution of ln(g 2 ) as a function of ln(L/15) for [to] = 2 years using Data 
set {2}. 

Figure 4. 

One example of a fit to the Data set {2} for [to] = 1 year and L = 1 day 
by a Pade Approximant (continuous curve). The chosen subset (see (19)) 
corresponds to 2000 points in the region Jomm = 420 and Jomax = 553. The Jo 
is computed to be 502 basis points. The empirical points are plotted versus 
the best Pade fit. In this case, the values of qi evaluated in their respective 
units are q 1 = 38.6 ± 1.3(%) _1 and q 2 = -308 ± 19(%)~ 2 . The obtained \ 2 is 
61 with 83 degrees of freedom. 

Figure 5. 

The parameter ln(g x ) as a function of Jo/600 for [m] = 1 year and L = 1 
day. These points have been obtained by minimizing the x 2 as a function of 
parameter q±. Each point corresponds to a subset of data containing close 
to two thousand points with average Jo value equal to Jo. The subsets are 
overlapping so that the points are not independent. The unit used for q 1 
is (%) _1 . For J it is the basis point. The straight line is a linear fit. Its 
equation is \n( qi ) = ln(/3i) - - //iln(l/15) with v x = 1.32 ± 0.02 

and ln(/3i) = 3.14 ±0.02. 

Figure 6. 
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The parameter ln(g 2 ) as a function of J /600 for [m] = 1 year and L = 1 
day. These points have been obtained using the relation defined in (16). Each 
point corresponds to a subset of data containing close to two thousand points 
with average J value equal to Iq. The subsets are overlapping so that the 
points are not independent. The error on ln(g 2 ) is computed by propagating 
the error on the estimate of the variance. The unit used for q 2 is (%)~ 2 . For 
Jo it is the basis point. The straight line is a linear fit. Its equation is ln(g 2 ) = 
ln(/3 2 )-i/ 2 (io/600)-// 2 ln(l/15) with v 2 = 2.19±0.02 and ln(/3 2 ) = 4.51±0.02. 

Figure 7. 

The parameter ln^) as a function of Jo/600 for [m] = 2, 3, 5, 7, 10, 30 
years and L = 1 day. These points have been obtained by minimizing the 
X 2 as a function of parameter q 1 . Each point corresponds to a subset of data 
containing close to two thousand points with average Jo value equal to Jo. The 
subsets are overlapping so that the points are not independent. The unit used 
for q 1 is For J it is the basis point. 

Figure 8. 

The parameter ln(g 2 ) as a function of Jo/600 for [m] = 2, 3, 5, 7, 10, 30 years 
and L = 1 day. These points have been obtained using the relation (16). Each 
point corresponds to a subset of data containing close to two thousand points 
with average J value equal to J . The subsets are overlapping so that the 
points are not independent. The unit used for g 2 is (%) -2 . For J it is the 
basis point. The error on ln(g 2 ) is computed by propagating the error on the 
estimate of the variance. 

Figure 9. 

Distribution of ln^), for [m] = 5 years as a function of J /600 for 3 different 
L values (1,10 and 30 days) where the L= 10 and 30 days distributions have 
been shifted on the ln^) axis by the correcting term (27). 

Figure 10. 

Distribution of ln(g 2 ), for [m] = 5 years as a function of J /600 for 3 different 
L values (1,10 and 30 days) where the L= 10 and 30 days distributions have 
been shifted on the ln(g 2 ) axis by the correcting term (27). 
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Table Caption 

Table 1. 

The values of the best fits to the parameters Hi (see (25)) using Data set {2}. 
These parameters are given as a function of the maturity [to] . The standard 
errors correspond to a unit deviation in the x 2 ■ 

Table 2. 

The values of the best fits to the parameters Hi (see (25)) using Data set {1}. 
These parameters are given as a function of the maturity [to] . The standard 
errors correspond to a unit deviation in the x 2 ■ 
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Figure (1) 

weight distribution for [m] = 1 year 




I ^P) 

Figure (2) 

ln(qj) versus ln(L) ([m] = 2 years) 
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Figure (3) 

ln(q 2 ) versus ln(L) ([m] = 2 years) 
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Figure (4) 

A fit for [m] = 1 year, L = 1 day and average L= 502 bp 
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Figure (5) 



ln(qj) versus L (L= 1 day & [m] = 1 year) 
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Figure (6) 



ln(q 2 ) versus I Q (L= 1 day & [m] = 1 year) 
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ln^j) versus L (L= 1 day & [m] = 2 years) 



Figure (7) 

s) & v (n 



(qj) versus I (L= 1 day & [m] = 3 years) 
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ln(qj) versus I (L= 1 day & [m] = 5 years) 



ln(q,) versus I„ (L= 1 day & [m] = 7 years) 
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ln(q[) versus I (L= 1 day & [m] = 10 years) 




ln(qj) versus I„ (L= 1 day & [m] = 30 years) 
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Figure (9) 

ln(qj) versus I ([m] = 5 years) 
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Figure (10) 

ln(q 2 ) versus L ([m] = 5 years) 
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Table (1) 

Data set {2} 

maturity (yrs) ^ (%~ r ) /i 2 (%~ 2 ) (%~ x ) 

1 0.657 ± 0.005 1.197 ± 0.005 1.82 ± 0.02 

2 0.627 ± 0.006 1.178 ± 0.007 1.88 ± 0.03 

3 0.653 ± 0.005 1.164 ± 0.006 1.78 ± 0.02 
5 0.638 ± 0.005 1.139 ± 0.005 1.79 ± 0.02 
7 0.626 ± 0.005 1.111 ± 0.005 1.78 ± 0.02 
10 0.636 ± 0.005 1.100 ± 0.005 1.73 ± 0.02 
30 0.578 ± 0.006 1.068 ± 0.006 1.85 ± 0.02 



Table (2) 

Data set {1} 

maturity (yrs) ^ (%~ x ) /i 2 (%~ 2 ) (%- x ) 

1 0.607 ± 0.006 1.209 ± 0.009 1.99 ± 0.03 

2 0.644 ± 0.004 1.200 ± 0.007 1.86 ± 0.03 

3 0.652 ± 0.004 1.172 ± 0.007 1.80 ± 0.02 
5 0.622 ± 0.006 1.139 ± 0.005 1.83 ± 0.02 
7 0.606 ± 0.006 1.109 ± 0.004 1.83 ± 0.02 
10 0.597 ± 0.005 1.093 ± 0.005 1.83 ± 0.02 
30 0.588 ± 0.005 1.074 ± 0.004 1.83 ± 0.02 
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